Linear and nonlinear information flow 
in spatially extended systems 
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Infinitesimal and finite amplitude error propagation in spatially extended systems are numerically 
and theoretically investigated. The information transport in these systems can be characterized in 
terms of the propagation velocity of perturbations Vp. A linear stability analysis is sufficient to 
capture all the relevant aspects associated to propagation of infinitesimal disturbances. In particular, 
this analysis gives the propagation velocity Vl of infinitesimal errors. If linear mechanisms prevail on 
the nonlinear ones Vp = Vl- On the contrary, if nonlinear effects are predominant finite amplitude 
disturbances can eventually propagate faster than infinitesimal ones (i.e. Vp > Vl)- The finite 
size Lyapunov exponent can be successfully employed to discriminate the linear or nonlinear origin 
of information flow. A generalization of finite size Lyapunov exponent to a comoving reference 
frame allows to state a marginal stability criterion able to provide Vp both in the linear and in the 
nonlinear case. Strong analogies are found between information spreading and propagation of fronts 
connecting steady states in reaction-diffusion systems. The analysis of the common characteristics 
of these two phenomena leads to a better understanding of the role played by linear and nonlinear 
mechanisms for the flow of information in spatially extended systems. 
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I. INTRODUCTION 

It is well recognized that chaotic dynamics generates 
a flow of information in bit space: due to the sensitive 
dependence on initial conditions one has an information 
flow from "insignificant" digits towards "significant" ones 
0. In spatially distributed systems, due to the spatial 
coupling, one has an information flow both in bit space 
and in real space. The flow in bit space is typically char- 
acterized in terms of the maximal growth A rate of in- 
finitesimal disturbances (i.e. of the maximal Lyapunov 
exponent) , while the spatial information flow can be mea- 
sured in terms of the maximal velocity of disturbance 
propagation Vp [Q-^. 

The evolution of a typical inflnitesimal disturbance in 
low dimensional systems is fully determined once the 
maximal Lyapunov exponent is known. The situation 
is more complicated in spatio-temporal chaotic systems, 
where infinitesimal perturbations can evolve both in time 
and in space. In this case a complete description of the 
dynamics in the tangent space requires the introduction 
of other indicators, e.g.: the comoving Lyapunov expo- 
nents and the spatial and the specific Lyapunov spec- 
tra §. 

Nevertheless, the complete knowledge of these Lya- 
punov spectra is not sufficient to fully characterize the 
irregular behaviors emerging in dynamical systems, this 
is particularly true when the evolution of finite perturba- 
tions is concerned. Indeed, finite disturbances, which are 
not confined in the tangent space, but are governed by the 
complete nonlinear dynamics, play a fundamental role in 
the erratic behaviors observed in some high dimensional 
system |^ A rather intriguing phenomenon, termed 



stable chaos, has been reported in j^]: the authors ob- 
served that even a linearly stable system (i.e. with A < 0) 
can display an erratic behavior with > 0. 

The first attempts to describe nonlinear perturbation 
evolution have been reported in p^ . However, in these 
studies the analysis was limited to the temporal growth 
rate associated with second order derivatives of one di- 
mensional maps. A considerable improvement along this 
direction has been recently achieved with the introduc- 
tion of the finite size Lyapunov exponent (FSLE) : a 
generalization of the maximal Lyapunov exponent able 
to describe also finite amplitude perturbation evolution. 
In particular, the FSLE has been already demonstrated 
useful in investigating high dimensional systems p^ ]. 

The aim of this paper is to fully characterize the in- 
finitesimal and finite amplitude perturbation evolution in 
spatio-temporal chaotic systems. Coupled map lattices 
(CML's) [|l5[ are employed to mimic spatially extended 
chaotic systems. The FSLE is successfully applied to 
discriminate the linear or nonlinear origin of informa- 
tion propagation in CML's. Moreover, a generalization 
of the FSLE to comoving reference frame (finite size co- 
moving Lyapunov exponent) allows to state a marginal 
stability criterion able to predict Vp in both cases: linear 
or nonlinear propagation. A parallel with front prop- 
agation in reaction-diffusion p6|jT7| ] (non-chaotic) sys- 
tems is worked out. The analogies between the two 
phenomena authorize to draw a correspondence between 
"pulled" ("pushed") fronts and linear (nonlinear) infor- 
mation spreading. 

The paper is organized as follows. In Sect. II the 
FSLE is introduced and applied to low dimensional sys- 
tems (i.e. to single chaotic maps). Sect. Ill is devoted 
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to the description and comparison of hnear and nonhn- 
ear disturbance propagation observed in different CML 
models. The finite size comoving Lyapunov exponent 
is introduced in Sect. IV and employed to introduce a 
generalized marginal stability criterion for the determi- 
nation of V^. A discussion on information propagation 
in non chaotic systems conclude Sect. IV. The analogies 
between disturbance propagation in chaotic systems and 
front propagation connecting steady states are analyzed 
in Sect. V. The Appendix is devoted to the estimation of 
finite time corrections for the computation of the FSLE 
in extended systems. Finally, some conclusive remarks 
are reported in Sect. VI. 



II. FINITE SIZE LYAPUNOV EXPONENT: 
LOW DIMENSIONAL MODELS 

Let us introduce the FSLE by considering the dynam- 
ical evolution of the state variable x = x(t) ruled by 

m = f [x(i)] , 

where f represents a chaotic flow in the phase space. 
In order to evaluate the growth rate of non infinitesi- 
mal perturbations one can proceed as follows: a refer- 
ence, x(t), and a perturbed, x (t), trajectories are con- 
sidered. The two orbits are initially placed at a dis- 
tance S{0) = Srnin, with Smiri ^ 1, assuming a cer- 
tain norm 6{t) = ||x {t) — x(i)||. In order to ensure 
that the perturbed orbit relaxes on the attractor a first 
scratch integration is performed for both the orbits un- 
til their distance has grown from 6m,in to Sq (where 
1 >> So >> Smi7i)- This transient ensures also the ahgn- 
ment of the perturbation along the direction of maximal 
expansion. Then the two trajectories are let to evolve 
and the growth of their distance 5{t) through different 
pre-assigned thresholds {5n = Sor"-, with n = 0,. . . ,N 
and typically 1 < r < 2) is analyzed. 

After the first threshold, Sq, is attained the times 
t((5„, r) required for 5{t) to grow from (5„ up to Sn+i are 
registered. When the largest threshold 6n (which should 
be obviously chosen smaller than the attractor size) is 
reached, the perturbed trajectory is rescaled to the ini- 
tial distance Smin from the reference one. 

By repeating the above procedure J\f times, for each 
threshold (5„, one obtains a set of "doubling" times (this 
terminology is strictly speaking correct only if r = 2) 
{Ti{Sn, f )}i=i,...,A'' and one can define the average of any 
observable A — A{t) on this set of doubling times as: 

AT 



where Ai = A{Ti{Sn,r)). The average (•)e does not coin- 
cide with an usual time average (•)* along a considered 
trajectory in the phase space, since the doubling times 



typically depend on the considered point along the tra- 
jectory and on the threshold The two averages are 
linked (at least in the continuous case) via the following 
straightforward relationship jl^ 



{A{t))t = - / dtAit) 



(At) 



(1) 



where T = J2i=i,M Tii^n.r) and (r((5„,r))e = T/Af. 

A natural definition of the Finite Size Lyapunov Ex- 
ponent A((5„) is the following |]l^: 



A(<S„) 
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(2) 



The last equality stems from the relationship among the 
two averages reported in (|l|). 

In the limit of infinitesimal perturbation (5„ and of infi- 
nite T (or Af) the FSLE converges to the usual maximal 
Lyapunov exponent 



lim lim A (5,,) = A . 



(3) 



In practice, at small enough Sn, A((5„) displays a plateau 
~ A. Moreover, one can verify that A((5„) is independent 
of r, at least for not too large r fl^ j. 

In Eq. (^) continuous time has been assumed, but dis- 
crete time is the most natural choice when experimental 
data sets (typically sampled at fixed intervals) are consid- 
ered. In order to generalize the FSLE's definition to the 
case of discrete time dynamical systems, let us consider 
the following map 

x(t + l) = F(x(i)), 

where x is a continuous variable, and t assumes inte- 
ger values. In this case T((5„,r) = t has simply to be 
interpreted as the minimum "integer" time such that 
(5(t) > Sn+i, and, since now 5{T)/5n is a fluctuating 
quantity, the following definition is obtained 



A(<5„) 
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(T((5„,r))c 
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5{t) 
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A theoretical estimation of (^) is rarely possible, and in 
most cases, one can only rely on a numerical computa- 
tion of \{5n). However, in the following we will report 
two simple cases for which an approximate analytic ex- 
pression for the FSLE can be worked out. 
Let us first consider the tent map 



F{x) 
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where x € [0 : 1]. This is a one dimensional chaotic map, 
since A = In 2 is positive. 

Due to the simplicity of this map, one can estimate 
the expression (|^) analytically obtaining the following 
approximation 
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A(5)~ln2-5, 



(5) 



valid for not too large S values. The maximal Lyapunov 
exponent is correctly recovered in the limit i5 — *■ and 
the above expression reproduces quite well the numerical 
estimate of the FSLE (see Fig. |l| (a) ). An important 
point to stress is that for this map the finite amplitude 
perturbations grow with the same rate or slower than 
the infinitesimal ones. The contraction of perturbations 
at large scales is due to saturation effects related to the 
attractor size. A similar dependence of the FSLE on 
the considered scale is observed for the majority of the 
chaotic maps (logistic, cubic, etc.), as we have verified. 
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FIG. 1. X{5) versus S for the tent map (a) and the shift 
map with l3 — 2 (b). The continuous lines are the analytically 
computed FSLE and the boxes the numerically evaluated one. 
The two maps are displayed in the insets. 

One can wonder if there are systems for which, at 
variance with the behavior (^), the finite size correc- 
tions leads to an enhancement of the growth rate at large 
scales. As shown in the shift map F{x) = Px mod 1 
represents a good candidate. Also in this case it is possi- 
ble to obtain an analytical expression for the FSLE, when 
<5< [l/(r-t-/3)]. 
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which again correctly reproduces the numerical data (see 
Fig. |l| (b) ) and in the limit 5^0 reduces to the cor- 
responding maximal Lyapunov exponent A = ln/3. As 
expected, finite amplitude disturbances can grow faster 
than infinitesimal ones : 



\{5) > A(<5 ^ 0) = A ior 0<5<5' 



(7) 



where (5^"* indicates the threshold at which saturation 
effects set in. An even more interesting situation is 
represented by the circle map F{x) = a + x mod 1. 
This map is marginally stable (i.e. A = 0), but it is 
unstable at finite scales. Indeed, the FSLE is given 
by X{5) = 5/{1 — 5)\n[{l — 5)/5], which is positive for 
< (5 < 1/2. Therefore at small, but finite, perturba- 
tions a positive growth rate is observed in spite of the 
(marginal) stability against infinitesimal perturbations. 
As a consequence, the circle map can exhibit behaviors 
that can be hardly distinguished from chaos under the 
influence of noise, since small perturbations may be occa- 
sionally driven into the nonlinear (unstable) regime and 
therefore amplified. Of course, the role of noise can be 
played by coupling with other maps, e.g. it has been 
found that coupled circle maps display behavior resem- 
bling (for some aspects) that of a chaotic system [ pl] |. 
This phenomenon becomes even more striking in certain 
coupled stable maps where, even if the maximal Lya- 
punov exponent is negative one can have a strong 
sensi tivity to non infinitesimal perturbations |l9[| (see 
Sect. IV B for a detailed discussion). 

The two maps here examined for which (Q) holds have 
a common characteristic: they are discontinuous. How- 
ever, in order to observe similar strong nonlinear effects, 
it is sufficient to consider a continuous map with high, 
but finite, first-derivative values ||l^. In this respect 
a simple example, that will be examined more in detail 
in Sect. IV, is represented by the map: 



hx 

F{x) = { l-c{l-q){x-l/b) 
q + dix-i^) 



0<x <l/b 
l/b<x<!^ 



b+c 
be 



< X< 



be 
1 



(8) 



with b = 2.7, d = 0.1, q = 0.07 and c = 500. For c ^ oo 
the map (||) reduces to the one studied in ||^. For the 
map (|^) the FSLE dependence on 5 is similar to that 
observed for the shift map. 



III. INFORMATION SPREADING IN SPATIALLY 
DISTRIBUTED SYSTEMS 



In this section we will examine the mechanisms behind 
the information flow in spatially distributed systems. In 
particular, the influence of linear and nonlinear effects 
on information (error) spreading will be analyzed. As a 
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prototype of spatially distributed system Coupled Map 
Lattices (CML's) |g5| are considered : 



Xilt) = {l-e)x,{t) + + .T,+i(<)) 



(9) 



where t and i are the discrete temporal and spatial in- 
dices, L is the lattice size (i — —L/2, . . . , L/2), Xi{t) the 
state variable, and £ S [0 : 1] measures the strength of 
the diffusive coupling. F{x) is a nonlinear map of the 
interval ruling the local dynamics. 

In order to understand how the information spreads 
along the chain, let us consider two replicas of the same 
system, x(i) = {xi{t)} and x'(i) = that initially 

differ only in a single site of the lattice (e.g. i = 0) of a 
quantity do, i.e. 



(10) 



The velocity Vp can be directly measured by detect- 
ing the leftmost, i;(<), and the rightmost, ir{t), sites for 
which at time t the perturbation Axi {t) exceeds a preas- 
signed threshold. The definition of Vp is the following 



Vr, = lim 



lim 



ir{t) - kit) 
2t 



(12) 



where the limit L — > oo has to be taken first to avoid 
boundaries effects. The velocity (|l3) does not depend on 
the chosen threshold values P,^,|lO|. 

Since the dynamics of the difference field ( ^t|) is not 
confined in the tangent space, non linearities can play a 
crucial role in the information propagation. Indeed, we 
will see that the evolution of the disturbances strongly 
depend on the considered map F{x) and in particular on 
the shape of A(i5). In the next subsection propagation in 
CML's with local chaotic maps for which X{6) < X V(5 is 
discussed. Local maps for which the condition (^) holds 
will be the subject of Sect. IIIB. 



where 5^ q is the Kronecker's delta. In a chaotic system 
the perturbation will typically grow locally and spread 
along the chain. These phenomena can be studied by 
considering the difference field 



Ax,{t) = \x'^{t) ~ x,{t)\ = \F{x.^{t - 1)) - F{S:,{t 



l))l- 
(11) 



It has to be stressed that the full nonlinear dynamics 
contributes to the evolution of Axi{t). 

The spreading of this initially localized disturbance can 
be characterized in terms of the velocity of information 
propagation Vp ||,^. As shown in Fig. ||, Axi{t) can 
grow only within a light-cone, determined by Vp. For ve- 
locities higher than Vp the disturbance is instead damped. 
This individuates a sort of predictability "horizon" in 
space-time: i.e. an interface separating the perturbed 
from the unperturbed region. 
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FIG. 2. Evolution of Axi{t), for a chain of coupled tent 
map lattices with a coupling e = 2/3. The initial perturba- 
tion is taken as in ( |lo| ) with do ~ 10~*. 



A. Linear Mechanisms 

Since in this subsection we consider CML's for which 
the local instabilities are essentially dominated by the 
behavior of infinitesimal perturbations, most of the fea- 
tures can be understood by limiting the analysis to the 
tangent space. 

The evolution in tangent space is obtained by lineariz- 
ing Eq. (I), i.e. 



5x^it + l) = F'{£,{t))[Sx,{t) + 
^{Sx^+i{t) ~ 2Sx,{t) + Sx,^i{t))] 



(13) 



where F' is the first derivative of a one dimensional 
chaotic map. Let us again consider as initial condition 
for the evolution of (^3|) a localized perturbation as ( p^ ) 
with do infinitesimal. The spatio-temporal dynamics of 
the tangent vector {dxi{t)} is determined by the interac- 
tion and competition of two different mechanisms present 
in Eq. (E^: the chaotic instability and the spatial diffu- 



As a first approximation, the effects of the two mech- 
anisms can be treated as independent. The chaotic in- 
stability leads to an average exponential growth of the 
infinitesimal disturbance: |5xo(t)| ~ doexp[Xt\. On the 
other hand, the spatial diffusion, due to the coupling, ap- 
proximately leads to a spatial Gaussian spreading of the 
disturbance : \6xi{t)\ « \6xo {t)\/\/ AnDt exp{-i'^ /iDt) 
where D = e/2. Combining these two effects one obtains 



\Sxi{t)\ « do 



exp 



2et / 



(14) 



Since the chaotic nature of the phenomenon will typically 
induces fiuctuations, Eq. (|lj) can only describe the aver- 
age shape of the disturbance. Moreover, Eq. ( p^ holds 
only when the perturbation is infinitesimal, since when 
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the disturbance reaches finite values a saturation mech- 
anism (due to the nonhnearities) sets in preventing the 
divergence of |i5a;i(t)|. 

To verify the vahdity of (p^, we studied the evolu- 
tion of localized perturbations of a homogeneous spatio- 
temporal chaotic state, in particular coupled logistic and 
tent maps have been considered in the regime of "fully de- 
veloped turbulence" j2|. Firstly, the system is randomly 
initialized and let to relax for a relatively long transient. 
At this stage two replicas of the same system are gener- 
ated and to one of the two a localized perturbation (like 
(p^)) is added. The evolution of the difference field (|ll|) is 
then monitored at successive times. In order to wash out 
the fluctuations, the shape of the disturbance is obtained 
averaging over many distinct realizations. 
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FIG. 3. Average evolution of perturbations for a CML of 
logistic maps {f{x) — 4a;(l — x)) for (a) e — 1/3 and (b) 
e = 1/10. {Axi{t)) is reported as a function of in a lin-log 
scale at different times (from bottom to top t — 10, 20, 30, 40). 
Deviations from a straight line correspond to deviation from 
the Gaussian shape. {Axi{t)) is obtained as an average over 
10^ realizations, for each one Aa:;i(0) has been chosen as ( [ic| ) 
with do = 10~^. For comparison the prediction (jl^) is also 
reported (dashed lines). 

As one can see from Fig. || Eq. (|lj) is fairly well veri- 
fied for large enough coupling while it fails at small e . 
These discrepancies are due to the finite spatial resolu- 
tion (that in CML's is always fixed to 1), since for small 
diffusivity constant the discretization of the Laplacian 
becomes inappropriate. The expression (^) for distur- 
bance evolution has been already proposed in Rcf. ||2|] 
for CML's in two dimensions. Deviations from (14) have 
been observed also in but attributed to anomalous 
diffusive behaviors. It has to be remarked that expres- 
sion (^4|) is valid only at short times, since asymptotically 
(t — > oo) the infinitesimal leading edge of the propagating 
front |(5xi(i)| assumes an exponential profile pd| |. 

For what concerns the propagation velocity, an esti- 
mation of Vp can be obtained for infinitesimal perturba- 
tions by the evaluation of the so-called maximal comov- 
ing Lyapunov exponents A{v) The time evolution of 



an intially localized (infinitesimal) disturbance ( [To| ) in a 
reference frame moving with velocity v can be expressed 



(15) 



by following the perturbation along the world line i — vt 
one can easily measures the corresponding comoving Lya- 
punov exponent A{v) (for more details see [||J2^]). The 
information propagation velocity is the maximal velocity 
for which a disturbance still propagates without being 
damped. Therefore it can be defined through the follow- 
ing marginal stability criterion : 



A(14) = 0, 



(16) 



where the velocity has been now indicated with Vl in or- 
der to stress that it has been obtained via a linear anal- 
ysis. For the maps considered in this section the identity 
Vp — Vl is always fulfilled. 




FIG. 4. Comparison between the directly measured prop- 
agation velocities Vp = Vl (circles) and the prediction ( |lq ) 
(boxes) for a CML of logistic maps with (a) a — 3.9 and (b) 
a = 4, and of tent maps with (c) a — 2 (d) a — 1.8. Lattices 
of 4 ■ 10* maps has been used. 

As shown in 1^ , in a closed system with symmetric cou- 
pling A = A(v) has typically a concave shape, with the 
maximum located at v = (in particular A = A(w — 0)). 
An approximate expression can be obtained for A(u), by 
substituting i — vt in (|lj) and by comparing it with ([l5| ) 



A{v) ^ X~v^/2e . (17) 

This parabolic expression for A(v) suffers of the same 
limits mentioned for the Gaussian approximation ( p^ ) 
for the disturbance evolution. Anyway, from (|l^) an an- 
alytical prediction can be obtained for Vl '■ 



Va = 



(18) 



which, as shown in Fig. ^, is indeed very good apart from 
some deviations for e ss and e 1. In Sect. M we will 
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rederive (|18|) by assuming that the chaotic perturbation 
behaves as a front connecting a stable and an unstable 
(metastable) fixed point in a non-chaotic reaction diffu- 
sion system. 

Let us briefly recall that another method (not suffer- 
ing for boundary problems) to determine the comoving 
Lyapunov exponent has been introduced in [ p2[ . The 
method relies on the computation of specific Lyapunov 
exponents A(/Lt) associated to an exponentially decaying 
perturbation (with spatial decay rate /z) . In other words 
one assumes that the spatio-temporal evolution of an in- 
finitesimal disturbance can be written as 



(19) 



Since the asymptotic leading edge of the front separat- 
ing perturbed from unperturbed part of the chain has an 
exponential shape, the above assumption ( p9[ ) is appro- 
priate to study its evolution. 

It is straightforward to show that the comoving Lya- 
punov exponents are related to the specific ones via a 



Legendre transform |22|, all the data concerning comov- 
ing exponents reported in this paper have been obtained 
with such a method. Moreover, a further results con- 
cerns the linear velocity Vl, it can be shown that 
its value corresponds to the minimal propagation veloc- 
ity VUi) = X{fi)/fi associated with perturbations of the 
form (119), i.e. 



Vl = mm — — 



(20) 



where /zl and Al = Mp-l) represent the spatial decay 
rate and the temporal growth rate of the leading edge, 
respectively. The expression ( |20| ) for the linear velocity 
is identical to the one derived for propagation of fronts 
connecting a stable to an unstable steady state |23] . 



B. Non Linear Mechanisms 

In this section we investigate the case of coupled maps 
for which X{6) > A(0) in some interval of 5. As noticed in 
Sect. H, this behavior can be observed in chaotic (abso- 
lutely unstable) maps, as well as in stable and marginally 
stable maps. Let us first analyze chaotic maps, non- 



chaotic ones will be discussed in Sect. IV B. 

For these systems it is possible to have Vp > Vl, this 
means that the disturbance can still propagate also in 
the velocity range [Vl , Vp] , even if the corresponding co- 
moving Lyapunov exponents are negative. Therefore, the 
linear marginal stability criterion ( |lq ) does not hold any- 
more. We want to stress that the condition (^ is neces- 
sary, but not sufficient to ensure that Vp > Vl, since all 
the details of the coupled model play also an important 
role. 

In Fig. ^ the spatio-temporal evolution of an initially 
localized disturbance of a chain of coupled shift maps is 



reported. As shown in |24j, when the coupling £ < 1/2 
the maximal Lyapunov exponent for such model coin- 
cides with that of the single map (namely, A = ln(/3)) 
and if /3 > 1 the system is chaotic. Initially the pertur- 
bation, that is still infinitesimal, spreads with the linear 
velocity Vl above defined. At later time it begins to 
propagate faster with a velocity Vp > Vl- Comparing 
Fig. I with Fig. |, one can see that the second stage of 
the propagation sets in when the bulk of the perturbed 
region reaches sufficiently high values. As a matter of 
fact the initial stage of propagation disappears if we ini- 
tialize the two replicas with a disturbance of amplitude 
0(1). From these facts it is evident that the origin of 
the information propagation characterized hy Vp > Vl 
should be due to the strong nonlinear effects present in 
this type of CML. 
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FIG. 5. The same of Fig. 
maps with l3 — 1.03, e — 1/3. 



for a lattice of coupled shift 



The behavior at long times can be understood by con- 
sidering the dependence of X{d) on the disturbance am- 
plitude as shown in Fig. actually the figure refers 
to the single map with (3 = 2, but the shape of X{6) is 
qualitatively the same also for the coupled system and 
for other /3-values. 

Until the perturbation is infinitesimal X{S) ~ A and 
the linear analysis applies, when the disturbance becomes 
bigger than a certain amplitude, 5^^ , the growth will be 
faster, since now X{6) > X. As it can be clearly seen in 
Fig. 1^, the perturbation is well reproduced by the linear 
approximation (|lj) until the amplitude of the perturba- 
tion reaches a critical value S^^ ^ 0(10"'') above which 
the nonlinear effects set in. At this stage the nonlinear 
instabilities begin to push the front leading to an increase 
of its velocity and deforming the profile of the perturba- 
tion. This becomes exponential at much shorter times 
than in the linear situation discussed in previous sub- 
section. Moreover, when the propagation is dominated 
by nonlinear mechanisms the spatial decay rate (j^nl of 
the asymptotic leading edge will be greater of the linear 
expected value fj,L'- this result can be explained again in- 
voking the analogy with propagation of fronts connecting 
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steady states [Tl| ]. 
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FIG. 6. Evolution of the perturbation Axi{t) for a CML of 
shift maps with e = 1/3 and /3 — 1.04 at four different times, 
t = 250, 450j650, 1000. The sohd hnes are the expected Gaus- 
sian shape (hj) . The decay rate of the asymptotic exponential 
profile is /inl ~ 1.47, noticeably greater than the linear value 
liL = 0.42. 
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FIG. 7. Time evolution of Axo{t) (solid line) and Axi{t) 
(dashed line), for i = V2t with 1/2 = 1/3. In this last case 
the data for two different chain configurations are reported. 
The map used is the shift map with 13 = 1.1, e = 1/3 and 
L — 2 ■ 10^. For these parameters one has Vl = 0.250 and 
Vp = 0.342. Note that between t = 200 and t = 500 for Ax^{t) 
the numerical precision is reached. 

To better clarify the difference between the linear and 
non linear mechanisms we show the behavior of Axi(t) 
along the world line i = Vt for two different velocities: 
Vl = and Vl <V2 < Vp. From Fig. |, one observes 
for the zero-velocity situation an exponential increase of 
the perturbation with rate A until it eventually saturates. 
For the case corresponding to velocity V2 an initial ex- 
ponential decay with rate A(V2) is seen, followed at later 
times by a resurgence of the disturbance. The successive 
evolution of the perturbation is no more exponential and 



exhibits strong fluctuations. 

These features suggest that in order to generalize the 
criterion (|l6| ) to nonlinear driven information spreading 
the growth of finite amplitude perturbations in a moving 
reference frame should be analyzed. 



IV. FINITE SIZE LYAPUNOV EXPONENT: 
EXTENDED SYSTEMS 

In this section we introduce the Finite Size Comov- 
ing Lyapunov Exponent (FSCLE) that is a generaliza- 
tion of the FSLE to a moving reference frame. First of 
all let us define the FSLE for an extended system: in 
this case exactly the same definition given in Sect. ap- 
plies, apart from some ambiguities in the choice of the 
norm to employ for measuring the distance S(t) of the 
perturbed, x'(t) — {x^{t)}, from the unperturbed rephca, 
x{t) = {xi{t)}. A natural choice could be to perturb ran- 
domly x(i) and to look for the doubling times associated 
to the evolution of the distance 

L/2 

^W = ;^ E ; (21) 

i=-L/2 

an alternative choice consists in perturbing a single site 
of the chain, let's say i = 0, at time t = and to evaluate 
the "single site" norm 

Axo{t) = \x'o{t) - xo{t)\ . (22) 

We have verified that the two norms ( |2l| ) and (|2^) give 
equivalent results for what concerns the evaluation of 
.^J). In this paper we will limit to consider the norm 

In order to measure the FSCLE in a reference frame 
moving with velocity v, we have simply measured the 
difference ( |2^ ) along a world line i = vt; i.e. 

Ax,^+[^t] (t) = Wi,+[vt] (t) - ^ia+lvt] (t) I • (23) 

where [ • ] denotes the integer part and ic is introduced 
below. 

The FSCLE is then estimated as in Sect. |l|. Once a 
set of thresholds (5„ = r"(So, with n = 0, . . . , A^, is chosen 
and the perturbation is initialized as Axi(O) — 5mm Sift 
with Smin ^ (5o. A preliminary transient evolution is per- 
formed in order to allow to the perturbed orbit to relax 
on the attractor. At the end of this short transient the 
position ic where the perturbation reaches its maximum 
is detected, this point is taken as the vertex of the light 
cone from which all the considered world lines i = vt 
depart. 

The FSCLE is then defined as 
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where the dependence on the velocity v derives from the 
employed norm . Note that we have used Eq. due 
to the discreteness of the temporal evolution. 

In the limit of very small perturbation the FSCLE re- 
duces to the comoving Lyapunov exponent 

lim A(J,w) = A(u), (25) 

and, ior V — one has the FSLE X{S). 

Actuallythere are finite time effects which prevents 
the limit (|2^) to be correctly attained. This is related to 
the fact that the FSCLE's can be obtained only via finite 
time measurements. In the Appendix we show how one 
can include such finite time corrections. 

In the next subsections we will give evidences that the 
marginal linear stability criterion ( p^ can be generalized 
with the aid of the FSCLE in the following way: 

max{ A((5, v)} = ioTv>Vp, (26) 

s 

where Vp can be either Vl, if the information propaga- 
tion is due to linear mechanisms, or greater than Vl, 
when nonlinear mechanisms prevail on the linear ones. 
For V > Vp one has A{6,v) — 0, since due to the defi- 
nition of the FSCLE negative growth rate appear to be 
0. 



In the interval Vl < v < Vp (see Fig. |(b)), A{5,v) 
is still positive and its maximal value decreases for in- 
creasing V. The FSCLE vanishes for v = Vp as Fig. ^ 
(a) shows. In this velocity range A{v) is always negative, 
therefore the instabilities observed in a reference frame 
moving with velocity higher than Vl have a fully non- 
linear origin. The residual fluctuations present in A{d, v) 
are essentially due to the limited statistics. Indeed the 
nonlinear growth as shown in Fig. |^ is extremely fluctu- 
ating. 
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A. Chaotic Systems 

Let us first consider chaotic systems for which Vp = Vl, 
in this case 

max{A((5, v)} = A{v) 
s 

and the generalized criterion (|2^) reduces to the linear 
one (p^). 

As already stressed in Sect. Ill B| , shift coupled maps 
represent a prototype of the class of chaotic models for 
which Vp can be eventually bigger than Vl- For this 
model the behavior of A{6,v) for various velocities is re- 
ported in Fig. ^. 

For t; < Vl, we observe that A{S,v) ~ A(t;) up to a 
certain value of the disturbance amplitude 6^^, above 
which the FSCLE increases and exhibits a clear peak at 
some higher 5- value. From Fig. || (a) it is clear that 
(which denotes the set in of the regime dominated by 
nonlinear mechanisms) decreases for increasing velocities 
and finally vanishes ed v — Vl- This behavior is reason- 
able since, as shown in Fig. |^, initially the disturbance 
evolves along the chain following the linear mechanism 
characterized by A{v), but as soon as in the central site 
5 > 5^^ the nonlinear mechanism begins to be active. 
Thus a nonlinear front is excited and this invades the 
linear region propagating with a velocity higher than Vl , 
therefore the higher is v the smaller is the scale at which 
nonlinear effects are observed. 
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FIG. 8. A((S, v) as a function 5 for various velocities. The 
data refer to the coupled shift maps with 13 = 1.1 and e = 1/3, 



suits for velocities v < Vl are reported, namely (from top to 
bottom) V — 0, V — 0.10 and v — 0.16. The straight lines 
indicate A(w). In (b) velocities in the range [Vl '■ Vp] axe re- 
ported, from the top v = 0.25, v = 0.286, v = 0.30, v = 0.33, 
V = 0.34. Chains of lengths from L = 10* up to L = lO'^ 
have been employed and the statistics is over 2 • 10^ doubling 
times. 

In Figure ^(a) the dependence of maxs{A{5,v)} on v 
for coupled shift maps is reported, as expected it vanishes 
exactly for v = Vp- For v < Vl & non-monotonous behav- 
ior of ma.xs{A{6,v)} is observable, this is probably due 
to the complex interplay of linear and nonlinear effects. 
For V > Vl, the behavior is smoother and a monotonous 
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decrease is observed. As discussed in Sect. ||, the dis- 
continuity present in the shift map is not necessary in 
order to observe nonhnear mechanisms prevaihng on hn- 
ear ones. 

As an example of continuous map exhibiting an in- 
formation propagation velocity Vp > Vl the map (||) is 
considered. This type of CML has been already studied 
in [l0|] : it has been observed that in a certain parameter 
range Vp can be finite even if the map is non chaotic. 
Moreover, also in the chaotic regime there is a window of 
parameters where Vp > Vl, with our parameters choice 
A « 0.182 > 0. 
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FIG. 9. maxg A{S, v) (dashed line witli points) versus v 
compared with A{v) (continuous line), (a) Results for coupled 
shift maps with the same parameters as Fig, m, the vertical 
lines indicates Vl ~ 0.250 and the measured (nJP propagation 
velocity Vp ~ 0.342. (b) Data for the coupled maps (H) with 
parameters b = 2.7, d = 0.1, q — 0.07 and c — 500 and with 
coupling e = 2/3, the vertical lines indicates Vl ~ 0.39 and 
Vp Ri 0.59. In this second case chains of length L — 2 ■ 10'* 
have been considered and the statistics is over 2-10*' doubling 
times. The reported values for max^ A{5, v) refer to an aver- 
age over five values of S around the peak position of A{S, v) . 



Also in the present case we observe an overall behavior 
of the FSCLE resembling that of the coupled shift maps. 
The main point that we want to remark is that in the 



limit V ^ Vp A{5,v) 0, as shown in Fig. ^ (b). At 
variance with the shift map (that is an everywhere ex- 
panding map) the considered map shows contracting and 
expanding intervals. Therefore the disturbances during 
their spatio-temporal evolution can be alternatively ex- 
panded or contracted. This leads to strong fluctuations 
in the FSCLE-values, that are difficult to remove. 

As a matter of fact we observed that even for veloci- 
ties slightly larger than Vp A(^, v) can be non zero. But 
the statistics of these anomalous fluctuations is extremely 
low, referring to the parameters reported in Fig. |^ (b) for 
V = 0.6 > Vp = 0.59 we observed an expansion (instead 
of the expected contraction) of disturbances in the 0.5% 
of the studied cases. For higher velocities A(5, v) is zero 
for any considered S. 



B. Non Chaotic Systems 

On the basis of the linear analysis discussed in 
Sect. Ill A the propagation of disturbances should be 
present only in chaotic systems. For non chaotic ones 
A < and from Eq. (|l^) one trivially obtains Vl — 0- 
On the other hand, in the class of systems for which 
A ((5) > A there are also stable and marginally stable sys- 
tems. Therefore, a propagation due solely to non linear 
terms can still be present [|9|-[ll|] . In this section we want 
to discuss the possible employ of the FSLE to character- 
ize these maps. 

Before entering into the description of the FSLE com- 
putation in such systems, it is of interest to recall an 
important phenomenon which appears in stable systems: 
the so-called "stable chaos" |^]. Stable systems asymp- 
totically evolve towards trivial attractors (i.e. fixed 
points, periodic or quasi-periodic orbits). However, in 
spatially extended systems it may happen that the time 
needed to reach the asymptotic state is very long: it 
has been found that in certain stable CML the transient 
time diverges exponentially with the number of elements 
of the chain Moreover, this transient is character- 
ized by a quasi-stationary behavior allowing for the in- 
vestigation of the properties of the model with statistical 
consistency. In Ref. ^ it has been shown that a chain 
of coupled maps of the type (||) , considered in their dis- 
continuous limit (i.e. for c — *■ oo), is non chaotic but 
still exhibits erratic behaviors. This is associated to a 
non-zero information spreading within the system. 

As far as the computation of the FSLE for these sys- 
tems is concerned, some remarks are worth to be done. 
The definition of the FSLE in term of error doubling 
times cannot be used in a straightforward manner to 
determine negative expansion rates. Another important 
point is that, at variance with the case of chaotic maps, 
finite perturbations should be now considered in order to 
observe an expansion. 

These two points impede a straightforward implemen- 
tation of our method to study these systems. Indeed, if 
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one starts with too small perturbations the propagation 
does not manifest, while if one initializes the system with 
a finite perturbation X{5) cannot be estimated with the 
required accuracy. Indeed, the only way to have an inde- 
pendence of A ((5) from the initial conditions is to initialize 
the system with infinitesimal perturbations and then to 
follow them until they become finite due to the dynamics 
of the system (see Ref. for a detailed discussion). 

The coupled maps for c — > oo have been recently 
analyzed by Letz & Kantz ||l9| in term of an indicator 
similar to the FSLE (i.e. able to quantify the growth rate 
of non infinitesimal perturbations). This indicator turns 
out to be negative for infinitesimal perturbations and be- 
comes positive for finite perturbations. This means that a 
finite perturbation of sufficient amplitude can propagate 
along the system due to nonlinear effects. This confirms 
what was previously observed in Ref. pl for marginally 
stable systems. 
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FIG. 10. I{S) as a function of the amplitude perturbation 
in a lin-Iog scale for the discontinuous coupled maps ^ , in the 
limit c — > 00, with £ = 1/3 and L = 35. The map is reported 
in the inset. The negative maximal Lyapunov A = —0.105 is 
recovered at small scales. For the computation details see the 
text. In the inset the single map is shown. 

In Fig. no we show the behavior of a quantity I{Sn) 
similar to X(5n) which has been obtained as follows. We 
considered two trajectories at an initial distance af- 
ter one time step evolution the distance S between the 
trajectory is measured. Then one of the two trajectory 
is rescaled at a distance Sn from the other, keeping the 
direction of the perturbation unchanged, and the proce- 
dure is repeated several times and for several values of 
Sn. Then we averaged ln((5/5„) over many different initial 
conditions obtaining I{S„). For Sn 0, this is nothing 
but the usual algorithm for computing the maximal Lya- 
punov exponent [ p5[ . 

As discussed in Ref. @ this method suffers of the 
problem that when S„ is finite one is not sure to sam- 
ple correctly the measure on the statistically stationary 
state. Indeed a finite perturbation will generically bring 
the trajectory out from the "attractor" . Nevertheless the 



result is in good agreement with the one obtained in |Tg| ] 
and confirms that at the origin of the perturbation prop- 
agation in this system there should be a mechanism very 
similar to the one discussed for the shift map. Further 
studies related to the "stable chaos" phenomenon have 
been recently performed [E6|. 



V. ANALOGIES WITH FRONT PROPAGATION 
IN REACTION DIFFUSIONS SYSTEMS 

The perturbation evolution in spatially distributed sys- 
tems can be described as the motion of an interface sep- 
arating perturbed from unperturbed regions. In this 
spirit, one can wonder if and to what extent it is possible 
to draw an analogy between the evolution of this kind of 
interface and the propagation of fronts connecting steady 
states in reaction diffusion systems. As already noticed 
in Ref. [pd| , the two phenomena display many similari- 
ties. In the following we will discuss the similarities and 
differences, in particular we will introduce a simple phe- 
nomenological model which can help us in highlighting 
the analogies. 

Let us start by recalling the basic features of fronts 
propagation in reaction diffusion systems with reference 
to the Fisher-Kolmogorov-Petrovsky-Piscounov (FKPP) 
equation 1 16 



dtO{z, t) = Ddl^e{z, t) + G{e{z, t)) , 



(27) 



where d{z, t) represents the concentration of a reactants 
which diffuses and reacts, the chemical kinetics is given 
by G{9). Typically the function GiO) e Ci[0,l] (with 
G(0) = G(l) = 0) exhibits one stable (6* = 1) and 
one unstable [9 = 0) fixed point. Once the system is 
prepared on the stable state {9{z) = OV2;), an initial 
(sufficiently steep) perturbation (e.g. a step function 
9{z, t — Q) — 9(z — zo) ) will give rise to a smooth front 
moving with a velocity Vp, that will connect the unstable 
and the stable fixed points: as a result the stable state 
will invade the unstable one. This equation admits many 
different traveling solutions that are typically character- 
ized by their propagation velocities, however for a suffi- 
ciently steep initial perturbation of the unstable state the 
selected front is unique and its velocity Vp is bounded in 
the interval [I4im, Vmax]^ where 



and 



Vmax D sup s 

V o<e<i I 



IG{9)\ 



(28) 



(29) 



If G{9) is concave supQ^g^i{G{9) /9} is attained at 9 = 
(i.e. supQ^g^i{G{9)/9} — G (0)) and the selected veloc- 
ity is always the minimum one: the front is "pulled" by 
the growth of infinitesimal perturbations of the unstable 
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state. Otherwise, if G{9) is convex one can observe a 
velocity Vp > Vmm- the front is now "pushed" by the 
growth of finite ampUtude perturbations. These results 
are known as the Aronson & Weinberger theorem ||2^] 
(for more details see Ref. [0). The velocity Vmin can 
be easily obtained by performing a linear analysis of (|27| ) 
and by employing a marginal stability criterion p^ . 

The subject of our analysis is the spreading of per- 
turbations in a chaotic media. In order to compare our 
case with the FKPP, one we should consider the time 
evolution of the difference of two chaotic trajectories 
{lS.Xi{tY\i=i^L- However, the nature of the two phases 
separated by the front is now different from the FKPP 
case. The interface separates an unstable (Ax^ = 0) 
state from a "statistically stable" one. With the term 
"statistically stable" we mean that behind the front, in 
the bulk of the perturbed region, /S.Xi{t) does not con- 
verge to a stable fixed point (as for the FKPP) but it 
fluctuates in a stationary way around an average value. 
This suggests that a good model for reproducing this 
dynamical evolution would be a FKPP with a stochas- 
tic kinetics However, many similarities with FKPP 
can be established by neglecting the chaotic fluctuations 
and considering the average shape of the front [|ll|. In 
practice, this can be done by averaging the perturbation 
evolution over many different initial conditions. 

Once the chaotic fluctuations are neglected, one can 
express the average perturbation growth in any site of 
the chain via the following mean-field approximation: 



Ui{t + 1) = cxp[A(-iti(t))] Uiit) , 



(30) 



where i, t are the discrete space and time indices and 
Mj = (1 — e)ui -f e/2(ui+i -|- Ui-i). Here, Ui{t) indicates 
the "average" Axi{t) and X{u) the corresponding FSLE 
(at least for positive growth rates). In the limit u — s- 
A(m) X, while, for u ~ 0{1), X{u) should reflect the 
saturation effects associated to the nonlinear map. In 
the infinitesimal limit {u — >■ 0) one essentially recovers 
the model discussed in Ref. |29|] . 

By pas sing to continuous variables is possible to show 
that Eq. ( |30| ) can be reduced to Eq. ( p7| ) (at least at the 
leading order). In order to transform Eq. ( pO| ) in its con- 
tinuous version, let us introduce infinitesimal spatial and 
temporal resolutions dx and dt, and assume that the dif- 
fusive scaling holds, i.e. dx^ = dt. Limiting to the first 
order expansion in dt (second order in dx), one easily 
shows that the continuous counterpart of Eq. (pOh is 



dtu = X{u)u+-dl^u, 

which is nothing but (|2^) with 

G{u) 
X{u) = , 



(31) 



(32) 



and D = s/2. As a consequence of the identity (|32|), one 
should observe a pulled dynamics for the chaotic front 



if A > A(u) Vu > 0, and a pushed one could occur 
only if max„{A(M)} > A. This can be considered as a 
reformulation of the Aronson & Weinberger theorem p^ ] 
in the context of information propagation. The velocity 
bounds (Eq), (29) can be now identified with Va = V2eX, 
i.e. Eq. dig), the first one and with 



,/2e max{A(w)} 



(33) 



the second one. We stress again that the lower velocity 
bound is indeed represented by Vl and that it coincides 
with Va only for sufficiently strong diffusive coupling s. 

Another interesting point that we can investigate via a 
numerical simulation of the mean field effective equation 
( pO| ) is the dependence of the propagation properties on 
the specific shape of A(it). For instance, as previously 
conjectured, we expect that if a monotonous decreasing 
X{u) is considered one should observe linear propagation, 
only. Indeed, numerical integrations of the model (|3^), 
with the choice X{u) = A — Bu, being A and B positive 
constants, show that Vp = Va = V^eA, where A corre- 
sponds to the maximal Lyapunov of the effective model. 
A generalization of this simple model would require to 
consider A and B as fluctuating quantities generated by 
suitable stochastic processes. But while for the choice of 
A one can have some hint lE9|, this is not the case for B. 
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FIG. 11. Propagation velocities (symbols) for the model 
(|^ with \{u) given by (M) as a function of i? = ln(r-2), with 
e — 1/3, A — in 1.1 and C = A. The threshold values are 
fixed to (5^^ = 10"^ and (5'"* = 0.53 and a chain of 4 • lO" 
sites has been used. The curves refer to Va = V^eA (solid 
line) and Vs — \PleB (dashed line). Note that in the lin- 
ear case, B < A, there is a perfect agreement between the 
measured velocity and the prediction. 

Let us now consider the non linear propagation case, 
as we have previously seen a necessary condition in order 
to have > Vl is that max„{A(M)} > A for some finite 
u. For what concerns the actual value of the selected ve- 
locity, this will depend on the value of the diffusive cou- 
pling and on the specific shape of A(u). In the following 
we will examine the dependence of Vp on two quantities 
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that characterize X{u): the difference | max^ A(5) — A| and 
the scale S^^ at which the non Unear effects set in. 

As a first example let us consider X{u) as a step func- 
tion: 



Xiu) 




< u < 



(34) 



where A, B and C are positive quantities, with B > A, 
while i5^^ and 5^"* are amplitude thresholds. The param- 
eter A is nothing but the Lyapunov exponent, B mimics 
the non linear terms leading to an enhancement of the 
growth rate, and the last term mimics the damping of 
the perturbation due to the saturation effects. In Fig. |ll| 
it is reported the behavior of the propagation velocity for 
the model (^0|) with A(u) given by (Q) for various values 
of B, once A and C are fixed. From the figure is evident 
that if B < A then Vp = Va (i.e. we are in the linear 
regime), while as soon as _B > A an increase of Vp with 
respect to Va is observed. 
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FIG. 12. Propagation velocities Ip as a function of S for 
the model (|^), Q with e = 1/3, A = In(l.l), B = ln(1.3), 
C = A, 5^* = 0.53. 
while the boxes to the case 5^"* — Sq 
correspond to Va = V2eA and Vs — 



The crosses refer to the case 5""* = 5o°* i 
The two lines 



In the whole examined parameter range Vp is always 
reasonably approximated by 14, but smaller. By increas- 
ing B one observes an increase of the difference between 
the measured velocity and the linear prediction. These 
results confirm that the condition (^ is indeed a neces- 
sary condition in order to observe nonlinear propagation 
of information. 

We will now investigate the role of 6^^ in determining 
the propagation velocity. It is quite obvious that model- 
ing the dynamics via the Eqs. ( pO|) and (Q) the nonlinear 
effects will disappear in the limit Snl ^ (5*"* and this is 
indeed confirmed by the simulations (see Fig. |l^ . In or- 
der to examine a less obvious situation, a modification of 
the expression (53) is also considered. In this second case 
gsat _ gsat _|_ gNT assumcd, therefore by varying 5^^ 
the extension of the amplitude interval over which the 



nonlinear mechanism is active will not be modified (be- 
ing fixed to (Jg"*). But also in this second case Vp ^ Va 
for increasing S^^ , this indicates that the smaller are the 
perturbation amplitudes affected by the nonlinear mech- 
anisms the stronger will be the nonlinear effect on the 
velocity : Vp Vs for S^^ 0. 

Typically, in generic CML's by varying a control pa- 
rameter both the difference | max^ X{S) — X\ as well as 
(5^^ will change. Moreover, even the definition of 5^^ 
for a continuous X{S) is not obvious. In order to under- 



stand the validity of the mean-field approximation ( 34 ) in 
a more realistic case we consider the shift map. In par- 
ticular, for X{u) we employed the analytical expression 
(^) valid for the single uncoupled map. 




FIG. 13. Information propagation velocities for the shift 
map with e — 1/3: (boxes) linear velocities Vl and (crosses) 
directly measured nonlinear ones. The two lines correspond 
to Va ( [l^ ) (dotted line) and to the propagation velocity for 
the model ( |3o| ) with A(u) given by (solid lin e). The dashed 
curve with asterisks is Va = \Jle max* {A ((5)}. 

The FSLE for coupled shift maps is actually differ- 
ent from (^), but we are neglecting correlations among 
different sites and effects due to the specific measure as- 
sociated to the system. Nevertheless, a numerical inte- 
gration of Eq. ( ^0| ) equipped with (||) reproduces semi- 
quantitatively the features observed for a lattice of cou- 
pled shift maps (see Fig^rSl). In particular, also the sim- 
ple mean- field model (jsCpwith the choice (||) is able to 
forecast the observed transition from nonlinear to linear 
behavior for /? — > 2. 

Concluding this section we can safely affirm that ( |30| ) 
is a reasonable model to mimic the perturbation evolu- 
tion at a mean field level, neglecting the spatio-temporal 
fiuctuations and correlations. In the same fashion X(u) 
can be considered as an "effective" non linear kinetics for 
the perturbation evolution. 
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VI. FINAL REMARKS 
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In this paper information (error) propagation in ex- 
tended chaotic systems has been studied in details. In 
particular, we have analyzed the relevance of linear and 
nonlinear mechanisms for the propagation phenomena 
in spatio-temporal chaotic coupled map lattices. Linear 
stability analysis is not always able to fully character- 
ize disturbance propagation. This is particularly true 
for (marginally) stable systems with strong nonlineari- 
ties, where finite size perturbations are responsible for 
information spreading in the system. When the non- 
linear effects prevail on the linear ones the propagation 
velocity of information Vp can be higher than the linear 
velocity Vl- A necessary condition for the occurrence 
of information spreading induced by nonlinear mecha- 
nisms has been expressed in terms of the Finite Size 
Lyapunov Exponent (FSLE). We have also shown the 
existence of strong analogies between error propagation 
and front propagation in reaction-diffusion models. In 
particular, the above mentioned necessary condition is 
analogous to the Aronson & Weinberger theorem |27| for 
front connecting stable and unstable steady states. In 
the linear and nonlinear case, the propagation velocity 
Vp can be identified via an unique marginal stability cri- 
terion involving Finite Size Lyapunov Exponents defined 
in a moving reference frame. This result generalizes the 
corresponding linear criterion expressed in terms of the 
maximal comoving Lyapunov exponents for the identifi- 
cation of Vl . 

These results can be of some interest for the synchro- 
nization and the control of extended systems. It has 
been recently shown that the synchronization of coupled 
extended systems is strongly influenced by nonlinear ef- 
fects. In particular, the synchronization time exponen- 
tially diverges with the system size, even in non-chaotic 
situations, provided that > pG]. We believe that 
in these systems the appropriate indicator to character- 
ize such transition would be the "transverse" Lyapunov 
exponent ||3|], once extended to finite scales. As far as 
control schemes are concerned, since they rely mainly on 
linear analysis ]3^ , new nonlinear methods (e.g. based 
on the concept of FSLEs) should be introduced in order 
to control the erratic behaviors due to fully nonlinear 
mechanisms. 

A further aspect that should be addressed in future 
work concerns the extension of the applicability of the 
FSLE also to linearly stable systems: a candidate in this 
respect could be the indicator recently introduced in . 

Finally, it is reasonable to expect that the present anal- 
ysis is not limited to discrete models but that it can 
be applied to continuous extended systems described in 
terms of PDE's, e.g. to the complex Ginzburg-Landau 
equation. 
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APPENDIX A: 

In any numerical computation of the Lyapunov expo- 
nents one is forced to use a finite time approximation for 
an infinite time limit. Nevertheless, provided that the 
convergence to the asymptotic value is fast enough, this 
is not a dramatic problem. As a matter of fact, in low 
dimensional systems very fast convergence to the asymp- 
totic value is usually observed. However, this problem 
manifest more dramatically in high dimensional systems, 
where the time to align along the direction of maximal 
expansion could be very long p9|J33y34|] . 



In the present case this problem is complicated by the 
fact that the FSLE is intrinsically a finite time indica- 
tor. Indeed the time a perturbation takes to grow from 
a value S to rS is finite unless (5 — > 0. However, for the 
CML models here analyzed and for initially localized per- 
turbation (|o|) it is possible to evaluate the corrections to 
apply to the FSCLE, estimated at finite time, in order to 
recover, for sufficiently small 6, the expected limit A(u). 

These corrections allow for a faster convergence of the 
FSCLE to its asymptotic values. 

For the sake of simplicity we consider maps with con- 
stant slope, e.g. the shift map F{x) — (3 x modi, and 
the CML defined in Eq. (|). In this case (for e < 1/2) the 
maximal Lyapunov exponent of the CML coincides with 
that of the single map |Q; i.e. A = ln/3. We will limit 
to the case v = Q (i.e to the FSLE), since the extension 
to generic v is straightforward. As shown in |3^,^ the 
finite time evolution of an infinitesimal perturbation do 
initially localized in i = can be expressed at time T as 
the sum of the contributions M{m, T) associated to all 
the paths connecting the space-time point (i = 0, t = 0) 
to the point (i = 0, t = T), i.e. 



fei=o(T) 
do 



^0 



T 

m=0 



M{m,T) 



(Al) 



where m is the number of "diagonal links" connecting 
{i,t) to (i ± 1,< + 1) present in the path of length T. 
Each path [m, T] of length T with m diagonal links con- 
tributes to the above sum with a term 



2) 



(A2) 
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where N{m,T) is the multipUcity associated to each 
path [to, T] (for more details see [Q). As can see froni 
Eq. (Al) the finite time Lyapunov exponent At will be 
given by 



At = ^ In 



5xo{T) 



do 



= ln/3+-ln 



^ Af(TO,T) 



ni—0 



(A3) 



the first contribution is the asymptotic one, while the 
second one will vanish in the limit T — s- cx) and it is the 
finite-time correction to evaluate. This second contri- 
bution can be numerically estimated by considering the 
finite time evolution of the Lyapunov eigenvector {Wi(t)} 
associated to the maximal Lyapunov exponent, once it is 
initiahzed as Wi{0) — doSi^o- The evolution of {^^^(i)} 
in the tangent space is ruled by the following equation 

W,{t + 1) = /3 [(1 - £)W,{t) + |(W^.-i(<) + W,+i{t)) 

(A4) 

In order to evaluate the finite time corrections, one 



should iterate at the same time Eq. (A4) (with f3 fixed 
to one) and the two replicas required for computing X{S) 
(see Sect. Then the estimation of the FSLE should 
modified in the following way: 



1 



(T((5„,r))c 



In 



In 



\6xi{t + T{dn,r))\ 

\Sx,{t)\ 

W^{t + T{Sn,r))\ 



(A5) 



The case of maps with non constant slope is computation- 
ally much heavier. Since for each different path the local 
multiplier F (xi) will be different and they will depend 
on the particular trajectory under consideration p3| , p5| . 
As a matter of fact we have observed that if the multipli- 
ers are equally distributed among positive and negative 
values, the finite time corrections essentially cancel out. 
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